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Abstract 

In this paper, the axial vibration of cracked beams, the free flexural vibrations 
of nanobeams and plates based on Timoshenko beam theory and first-order 
^ ■ shear deformable plate theory, respectively, using Eringen's nonlocal elas- 

ticity theory is numerically studied. The field variable is approximated by 
Lagrange polynomials and non-uniform rational B-splines. The influence of 
the nonlocal parameter, the beam and the plate aspect ratio and the bound- 
ary conditions on the natural frequency is numerically studied. The influence 
of a crack on axial vibration is also studied. The results obtained from this 
study are found to be in good agreement with those reported in the literature. 
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1. Introduction 



The inherent assumption in local elasticity is that dimensions of engineer- 
ing structure are so much larger than the characteristic dimensions of the 
microstructure. Thus, the classical continuum theories lack the capability 
of representing the size effects since they do not include any internal length 
scale |i|. Consequently, these theories fail when the specimen size or the 
wavelength become comparable with the internal length scales of the mate- 
rial. Several modifications of the classical elasticity formulation have been 
proposed, such as the strain gradient theory J2|, modified coupled stress the- 
ory |3|, |4| and nonlocal elasticity theory [5|, [g, |7|. A common feature of these 
theories is that they include one or several intrinsic length scales. The pre- 
dictions of these theories reduce to those of local continuum theories where 
the specimen size is much larger than the internal length scale. The key idea 
of the nonlocal elastic approaches is that within a nonlocal elastic medium, 
the particles influence one another not simply by contact forces and heat dif- 
fusion but also by long range cohesive forces. In this way, the internal length 
scale, can be considered in the constitutive equations simply as a material pa- 
rameter, called nonlocal parameter. One approach to estimate these nonlocal 
material parameters is by matching the phonon dispersion relation computed 
by these theories with the lattice dynamics dispersion relation [8|, |9|, Il0| . 

It is seen from the literature that the amount of work carried on the ap- 
plication of nonlocal and/or gradient elasticity theor y to study the dynamic 
response of nanobeams is considerable ll|, |l2|, [l3|, |l4], |l5] . Eringen's nonlocal 
elasticity theory has been applied by several authors to study axial vibra- 
tions 16[ and free transverse vibrations of nanostructures [12|, |T7|, [14( . Reddy 
and Pang 18( derived governing equations of motion for Euler-Bernoulli 
beams and Timoshenko beams using the nonlocal differential relations of 
Eringen [5j and presented closed form solutions for beam static bending, vi- 
bration and buckling response with various boundary conditions. Recently, 



Reddy [19|, based on von Karman nonlinear strains and Eringen's nonlo- 



cal theory, derived governing equations of equilibrium for beams. Phadikar 



and Pradhan [14| developed a variational formulation for Euler-Bernoulli 



beam and Kirchoff plate and employed Hermite and Lagrange polynomials to 
study the response of beams and plates within nonlocal elasticity framework. 
Roque et al., [15] studied bending, buckling and vibration of Timoshenko 



nanobeams using collocation techniques. Civalek and Demi [20( applied non 



local Euler- Bernoulli beam theory to study static bending and vibration of 



microtubules using differential quadrature method. Danesh et al, [21 



em- 



ployed differential quadrature method to study axial vibration of nanorod. 
The nonlocal linear elasticity of Eringen has been applied to study free in- 



plane vibration and flexural vibration of plates |22l . |23| . |24| . Reddy |l9j and 
Aghabahaei and Reddy 25| extended classical, first and third order shear 
deformation theory of plates to account for nonlocal effects. Aghababaei and 
Reddy 25[ reformulated the third-order shear deformation theory using non- 
local linear theory of elasticity and presented analytical solutions of bending 
and vibration of plates. The free vibration of single-layered and multi-layered 
graphene sheets have been studied using first order shear deformation non- 



local plate theory |23l . |26 



Nanostructures, like any other structures, may develop a flaw during man- 
ufacturing. For example, thermally induced cracks are observed in ZnO 
nanorods during fabrication 27|.The presence of a crack affects the dynamic 
response of a structural member. This is because, the presence of the crack 
introduces local flexibility. Moreover, with diminishing external characteris- 
tic length, the nonlocal effects cannot be ignored. Recently, Loya et al., 
and Hasheminejad et al., 
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studied transverse vibrations of nanobeams 
based on Euler-Bernoulli theory with and without surface effects, respec- 
tively. Hsu et al., [30] studied axial vibration of cracked nanobeam using 
nonlocal elasticity theory. 

To author's knowledge, the existing numerical approaches are limitedjx) us- 

to 
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ing radial basis functions [15| and differential quadrature method 
study the response of nanostructures based on nonlocal elasticity theory. 
In this paper, we use Lagrange basis functions and non-uniform rational B 
splines to approximate the field variable and the response is studied within 
the framework of the finite element method. The influence of the nonlocal 
parameter, the boundary conditions and the beam/plate aspect ratio on the 
fundamental frequency is numerically studied. The free axial vibrations of 
cracked nanobeams is studied by applying the extended finite element method 
(XFEM), where the crack is modelled independent of the mesh. This is done 
by adding suitable 'enrichment' functions to the already existing finite ele- 
ment basis functions. The influence of the boundary conditions, the crack 
location, the crack severity and the nonlocal parameter on the longitudinal 
frequency is studied. 



The paper is organized as follows. A brief introduction to Eringen's nonlocal 
elasticity is given in the next section, followed by the governing equations 
of motion for nanorod, nanobeam based on Timoshenko beam theory and 
nanoplate based on the first-order shear deformable plate theory (FSDT). 
Section [3] presents an overview of the basis functions used in this study and 
the corresponding discretized form. Section H] presents results for the free 
axial vibrations of cracked beam and free flexural vibration of Timoshenko 
beams and first-order shear deformable plates based on nonlocal elasticity 
theory. The effect of nonlocal parameter, the beam/plate aspect ratio and the 
boundary conditions on the fundamental frequency is numerically studied, 
following by concluding remarks in the last section. 

2. Nonlocal elasticity theory 

Eringen J6|, |7|, by including the long range cohesive forces, proposed a nonlo- 
cal theory in which the stress at a point depends on the strains of an extended 
region around that point in the body. Thus, the nonlocal stress tensor a at 
a point x is expressed as: 

<x= /"«(|x , -x|,r)t(x') dx' (1) 

where t(x) is the classical, macroscopic stress tensor at a point x and the 
kernel function a(|x' — x|) is called the nonlocal modulus, also referred to as 
the attenuation or influence function. The nonlocality effects at a reference 
point x produced by the strains at x and x' are included in the constitutive 
law by this function. The kernel weights the effects of the surrounding stress 
states. From the structure of the nonlocal constitutive equation, given by 
Equation (JTJ), it can be seen that the nonlocal modulus has the dimension of 
(length) -3 p, |7|. Typically, the kernel is a function of the Euclidean distance 
between the points x and x'. The kernel in Equation ([TJ) has the following 
properties: 

• The kernel is maximum at x' = x and decays with |x' — x|. 

• The classical elasticity limit is included in the limit of vanishing internal 



characteristic length 

lima(|x'-x|) =(J(|x'-x|) (2) 

T— >0 

• For small internal characteristic length, the nonlocal theory approxi- 
mates the atomic lattice dynamics. 

• The kernel is a Green's function of a linear differential operator, L 

M|x'-x|)=5(|x'-x|) (3) 

By appropriate choice of the kernel function, Eringen |5| showed that the 
nonlocal constitute equation given in integral form (Equation (pQ)) can be 
represented in an equivalent differential form as: 



(1 - t 2 L 2 V 2 )<t = t (4) 

where r 2 = j^ = (^ff ) , e is a material constant and a and L are the internal 
and external characteristic lengths, respectively. 

2.1. Rod 



For nanorod, the constitutive relation is given by [21 



d 2 N _ du 
dx 2 dx 

where EA is the effective axial rigidity and u is the axial displacement. The 
equations of motion for the axial vibration of nanorod is given by: 

9 fr, A 9u\ r d 2 \ d 2 u 

m \ EA Tx) = { 1 - ^) p w (6) 

2.2. Timoshenko beam theory 



For Timoshenko beam, the constitutive relations are given by 19 



M - /x 



d 2 M 



Q-fi 



dx 2 



dx 2 
GAk 



dx 
dw\ 

dx J 



(7) 



where Q is the shear force, G is the shear modulus and k s is the shear 
correction factor. The displacement field based on the Timoshenko beam 
theory is given by: 



u(x, z, t) = u (x, t) + z(f>(x, t); w(x, z, t) = w (x, t). 



(8) 



where u and w are the axial and the transverse displacements of the point 
on the mid-plane (i.e. z — 0) of the beam. The nonzero strain of the 
Timoshenko beam theory is 



duo (ty 
dx dx' 



7a 



dw 
dx 



+ 



(9) 



The Timoshenko beam theory requires shear correction factors to compensate 
for the error due to the constant shear stress assumption. The governing 
equations of motion in the case of nonlocal elasticity is given by: 



d_ 

dx 



^-[EI^\ -(7 4,,, 
dx \ dx 



GAk, [6 + — 



dw\ 

dx J 

dw 



+ 



dx 
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d 2 w 
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d 2 4> 



\x 



d 4 w 

dx 2 dt 2 

dU 



V 



dx 2 dt* 



(10) 



2.3. Reissner-Mindlin plate 

Consider a rectangular Cartesian coordinate system (x, y, z) with xy— plane 
coinciding with the undeformed middle plane of the plate and the z- coordi- 
nate is taken positive downwards. The displacement field based on first order 
shear deformation plate theory (FSDT), also referred to as the Mindlin plate 
theory is given by: 



u{x,y,z,t) 
v(x,y,z,t) 

w(x,y,z } t) 



u (x,y,t) +z9 x {x,y,t) 

v (x,y,t) + z6 y (x,y,t) 
w a (x,y,t) 



(11) 







zs b 



(12) 



The midplane strains e p , bending strain £{,, shear strain e s in Equation (fT2l 
are written as 



u 



o.x 



J o,y 



^o.y ~t~ Uq,x 



> , e b = < 
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o,, 



x,x 



'v,y 

"x,y * "y,x 



> e. 




(13) 

where the subscript 'comma' represents the partial derivative with respect 
to the spatial coordinate succeeding it. The equations of motion is given by: 



dN rr dN, 
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(14) 



h/2 

where (hi, h 2 , I 22 ) = J p(l,z,z 2 ) dz. The nonlocal membrane stress re- 

-h/2 

sultants N and the bending stress resultants M can be related to the mem- 



brane strains, e p and bending strains e b through the following constitutive 



relations 19 



N, 



j-.r 



N = { N yy } (1 - /"V 2 ) = A e e p + B be e b 
N xy 

M xx 
M { M yy } (1 - /iV 2 ) = B be e p + B b e b (15) 

M xy 

where the matrices A e = Ay, Bfe e = By and D& = Dy; (i, j = 1, 2, 6) are the 
extensional, bending-extensional coupling and bending stiffness coefficients 
and are defined as 



rh/2 _ 
{A ij} B iv D {] }= / Q y {l, 2, z 2 } dz, i = 1,2,6 (16) 



Similarly, the transverse shear force Q = {Qxz,Qyz} is related to the trans- 
verse shear strains e s through the following equation 

Qij = E lj£s (17) 

where By = E = j_ h/2 Qv(Uj dz; (i,j = 4, 5) is the transverse shear stiffness 
coefficient, Vi,Vj is the transverse shear correction factors for non-uniform 
shear strain distribution through the plate thickness. The stiffness coeffi- 
cients Q t j are defined as 



Q -Q -^- Q ~^- Q -Q -0 

V44 = Qbb = V66 = 2(1 _i_ 7y\ ( ) 

The equations of motion in terms of midplane displacements and rotations 
are obtained by substituting the nonlocal constitutive relations (Equation 




(TT5|) ) into the equations of motion Equation (fT4"|). For detailed derivation, 
interested readers are referred to the literature |2i 



3. Basis functions and spatial discretization 

In this study, the finite element model has been developed using Lagrange 
polynomials, NURBS basis function. In this section, a brief overview of 
different basis functions used in this study is given. As it is not the scope 
of this study to review the different basis functions, interested readers are 
referred to the literature and references therein. 



Lagrange polynomials. The Lagrange polynomials [3l| for a two noded ele- 
ment is given by: 

tfi = i-f; n 2 = ^ (19) 

where P. is the length of the element. The Lagrange polynomials are used to 
model the nanorod. 

Field consistent Q4 and Q8 element. The plate element employed here is a 
C° continuous shear flexible field consistent element with five degrees of free- 
dom (u , v , w , B x , By) at four nodes in an 4-noded quadrilateral (Q4) element 
and eight nodes in a 8- noded quadrilateral element (Q8). If the interpola- 
tion functions for Q4/Q8 are used directly to interpolate the five variables 
(u , v , w , B x , By) in deriving the shear strains and membrane strains, the el- 
ement will lock and show oscillations in the shear and membrane stresses. 
The field consistency requires that the transverse shear strains and mem- 
brane strains must be interpolated in a consistent manner. Thus, the B x and 
By terms in the expressions for shear strain e s have to be consistent with the 
derivative of the field functions, w 0>x and w 0>y . This is achieved by using field 
redistributed substitute shape functions to interpolate those specific terms, 



which must be consistent as described in 32, 331. This element is free from 



locking and has good convergence properties. For complete description of 



the element, interested readers are referred to the literature [32|, |33j], where 
the element behaviour is discussed in great detail. Since the element is based 
on the field consistency approach, exact integration is applied for calculating 
various strain energy terms. 



Non-uniform Rational B-Splines. The key ingredients in the construction of 
NURBS basis functions are: the knot vector (a non decreasing sequence of 
parameter value, £j < £i+i,i = 0, 1, ■ ■ • ,m — 1), the control points, Pi, the 
degree of the curve p and the weight associated to a control point, w. The 



i th B-spline basis function of degree p, denoted by N iyP is defined as [34J, |35| : 

1 if &<£<& + i 



N »® - 1 else 

N i>p (0 = -L^L-N^it) + , Wl ~/ ^ + i„-i(0 (20) 

Si+P ~~ ?i Si+p+1 — «+l 

A p' ft degree NURBS curve is defined as follows: 

m 

C(0 = ^ (21) 

i=0 

where Pj are the control points and Wi are the associated weights. 

The displacement field variable is approximated by one of the above men- 
tioned basis functions. By multiplying the governing equations of motion by 
appropriate test functions and by employing Bubnov-Galerkin method, we 
get the following set of algebraic equations: 

MS + K6 = (22) 

where M is the consistent mass matrix, K is the stiffness matrix and 8 is the 
vector of the degree of freedom associated to the displacement field. After 
substituting the characteristic of the time function 8 = u 2 8, the following 
algebraic equation is obtained: 

(K-lo 2 M)8 = (23) 

where u is the natural frequency. The frequencies are obtained using the 
standard generalized eigenvalue algorithm. 



10 



4. Numerical examples 

In this section, we present the axial vibrations of cracked beam and natu- 
ral frequencies for Timoshenko beam and first-order shear deformable plate, 
based on Eringen's nonlocal elasticity theory. In all cases, we present the 
non dimensionalized free flexural frequencies as, unless specified otherwise: 

For axial vibration of beams. 



{} = ^ L \I^A (24) 

where L, p and EA denote the length of the beam, mass density and effective 
axial rigidity. 

For flexural vibration of Timoshenko beam. 



-2./M 



^ = "L^>- (25) 

where EI is the flexural rigidity and p is the mass density. 

For plates. 

n = whj^ (26) 

where u is the natural frequency. In all cases, the effect of the nonlocal 
parameter is studied by the frequency ratio, which is defined as: 

Frequency ratio = — — (27) 

The frequency ratio is a measure of the error made by neglecting the nonlocal 
effects. For example, if -^- = 1, the effect of nonlocal parameter is not 
significant and becomes important for any other value. 

4-1. Axial vibration of cracked beams 

Consider a slender nanorod with uniform cross section A along the length L. 
The crack at distance of C from the left end is simulated by an equivalent 

11 



linear spring. The analytical expression for the eigen frequencies are given 
by 0: 




KPy/l-nP 2 cos 



b0 



'1 - ^ 



sin 



P-bfi 



(28) 



Clamped- Clamped. 

2 sin | , P 



+K/3 v / l-/i/3 2 



COS 



/9 



v/W^ 



cos 



b-2bp 



aA 7 ^ 



(29) 



when if (= -ff-), & and ju(= (^f 5 ) ) are known, the dimensionless longitudinal 
frequency Q of a single cracked nanobeam with nonlocal effect is solved by 
using Equations (128]) - (1291) . For a nanobeam without a crack, i.e., K = 0, 



the various vibration modes are expressed as [la, |36|, (30 



Clamped- free. 



/?: 



NL 



(2ra 



!)! 



1 + (2n- l) 2 /i(f)' 



:, 71= 1,2,3, 



(30) 



Clamped- clamped. 



P. 



mx 



NL 



1 + n 2 //. (7r) 



n = 1,2,3, 



(31) 



By choosing /i = in Equations (1301) and (l3Tj) . the various vibration modes for 
local elasticity are obtained. Figure (JT|) shows the convergence of first three 
modes with mesh size with and without nonlocal effects. It can be seen that 
with decreasing mesh size, the analytical solution is approached. Based on 
a progressive mesh refinement, 100 elements is found to be adequate for the 
study. Tables [H- [2] gives a comparison of computed frequencies for clamped- 
free and clamped-clamped nanobeam with a single crack, respectively. It can 
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be seen that the numerical results from the present formulation are found 
to be in good agreement with the existing solutions. The effect of nonlocal 
parameter on the fundamental frequencies is shown in Figure ([2]) for a fixed 
nonlocal parameter e a =1 nm. In this case, the length of the beam is varied. 
It can be seen that as the length of the beam increases with respect to the 
characteristic internal length, the nonlocal effects decreases and for very large 
beams, the nonlocal effects are negligible. 

Table 1: Comparison of natural frequencies of the clamped- free beam with a single crack 
located at C/L = 0.2002, crack severity K = 0.1144 and nonlocal parameter e a/L = 0. 

Mode XFEM Ref. [3pJ Ref. [37] 

1 1.4228 1.4278 1.4278 

2 4.4429 4.5579 4.5576 

3 7.8559 7.8540 8.8540 

4 10.4289 10.4471 10.4486 



Table 2: Comparison of natural frequencies of the clamped-clamped beam with a single 
crack located at C/L = 0.25 for various internal length and crack parameter [30J. 



e a/L 


K = 


0.065 


K = 


: 0.35 


K 


= 2 




XFEM 


Ref. [30] 


XFEM 


Ref. [30] 


XFEM 


Ref. [30] 


0.2 
0.4 


2.6144 
1.9455 


2.6173 
1.9467 


2.4649 
1.9060 


2.4668 
1.9071 


2.1503 
1.7660 


2.1506 
1.7663 



Figure ([3]) depicts the influence of boundary conditions and nonlocal param- 
eter on the first three modes of the clamped-free and clamped-clamped beam 
with a single crack. The crack is located at C/L = 0.4 from the left end. It 
can be seen that the frequencies are higher for the beam with no crack and 
with increasing crack parameter K the frequency decreases. The presence of 
the crack introduces local flexibility and the effect of crack parameter on the 
frequencies is significant for lower values of the nonlocal parameter e a. This 



observation is consistent with the existing results [30]. The effect of crack 



location along the length of the beam and the boundary conditions on the 

13 



10' 



o 

gio 

3 
0) 



.E 10' 



> 



3 10 



10' 









V 


2 




— h — mode 1 

mode 2 

mode 3 


* 


3 
4 
5 









10" 



10" 
h 



(a) e a/L = 



10" 



10 



10" 







— h — mode 1 


O _ 2 

g 10 




mode 2 

mode 3 


d) 




S 10- 3 




o 




elative 

o 

1 








10 10 10 

h 

(b) e a/L = 0.1 
Figure 1: Convergence of the first three modes with mesh size. 
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Figure 2: Frequency ratio ((3l/ Pni) as a function of beam length with the nonlocal param- 
eter e a = lnm. It can be seen that as the length of the beam increases with respect to 
the characteristic internal length, the nonlocal effects decreases and for very large beams, 
the nonlocal effects arc negligible. -. r 
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Figure 3: The dimensionless frequency for different crack parameters for the first three 
vibration modes. The crack is located at a distance 0.4 from the left end. 
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frequency is shown in Figure (Jl|). In both cases, the natural frequency of the 
beam is influenced by the location of the crack. In the case of clamped-free 
boundary condition, the crack near the free end shows a stronger influence 
than the crack at the fixed end, while, in the case of clamped-clamped bound- 
ary condition, the crack in the middle of the beam has stronger influence. 
Due to the symmetric boundary conditions, the frequency distribution in 
the case of clamped-clamped boundary condition is observed. This again is 
consistent with the results available in the literature [30|. 

4-2. Free flexural vibration of beams 

In the next example, consider a nanobeam of length a = 10 and thickness 
h and with following material properties: Young's modulus, E = 30 x 10 6 , 
Poisson's ratio v = 0.3 and mass density p — 1. The value of the shear 
correction factor k s is taken to be 5/6. 



Table 3: Comparison of natural frequencies of simply supported beams based on Timo- 
shcnko beam theory: °[12], ^1501 global collocation points [15], b NURBS basis functions 
with p = 3 and 82 elements, *Lagrange polynomials (100 elements). 
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FEM* 
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8.0750 


8.0423 


8.0748 


8.0706 







9.8381 


9.8059 


9.8281 
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9.3858 


9.3551 
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9.3452 
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8.0503 


8.0239 


8.0421 


8.0154 







9.7454 


9.7792 


9.7075 


9.5886 


10 
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9.2973 


9.3294 


9.2612 


9.1460 
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7.9744 


8.0014 


7.9434 


7.8445 



The numerical results for free vibration of Timoshenko beam are given in 
Table [3j It can be seen that the numerical results from the present study 
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Figure 4: The dimensionless frequency of mode 1 of (a) clamped-free beam and (b) 
clamped-clamped nanobcam. 
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are in good agreement with the existing solutions 12|, [l5| . Figure §5§ shows 
the convergence in case of the NURBS basis function, for mode 1 frequency 
with decreasing mesh size for various order of the curve. It can be seen 
that increasing the order of the curve, yields more accurate results. With 
decreasing mesh size, all the results approach analytical solution. Note that, 
the accuracy on a coarse mesh can be improved by increasing the order of the 
curve. The effect of nonlocal parameter and beam aspect ratio is shown in 
Figures flS} and © , respectively. It can be observed that increasing the non- 
local parameter decreases the fundamental frequency for both Euler-Bernoulli 
beam and Timoshenko beam. The nonlocal parameter has greater influence 
on higher modes. The variation of frequency ratio VL^l/^l is shown in Fig- 
ure ^. For mode 1, the variation is almost linear, but for higher modes, the 
variation becomes nonlinear. It can be seen from Figure ([7|) that increasing 
the beam aspect ratio increases the frequencies, while increasing the nonlocal 
parameter decreases the frequency. This is because, increasing the aspect ra- 
tio, increases the flexibility of the structure, on the other hand, the nonlocal 
parameter decreases the flexibility. The effect of various boundary condition 
on the fundamental frequency is illustrated in Table HI It can be seen that 
the effect of nonlocal parameter is to decrease the frequency. The clamped- 
clamped boundary condition increases the fundamental frequency, while the 
clamped-free lowers the frequency compared to the simply supported case. 




(total number of dots) 

Figure 5: Convergence of the first mode with mesh size for various degree of the curve. 
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Figure 6: Frequency ratio (fiiVL/^z,) as a function of nonlocal parameter /i for a/h = 
100 and L = 10 for Timoshenko beam. It can be seen that as the nonlocal parameter 
increases, the nonlocal effects are prominent. For a fixed nonlocal parameter, the effect of 
nonlocal parameter on higher frequencies is greater than that for lower modes. 
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Figure 7: Mode 1 non-dimensionalized frequency (Qnl) as a function of beam aspect ratio 
L/h for different nonlocal parameters for Timoshenko beam. It can be seen that the effect 
of nonlocal parameter is to decrease the frequency, whereas, increasing the beam aspect 
ratio, increases the fundamental frequency 
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Table 4: Non-dimensional frequency of Timoshcnko beam under various boundary condi- 
tions. 



a/h \x NURBS, p = 3 



SS CC CF 

9.8680 22.3892 3.5178 

1 9.4144 21.1228 3.4385 
100 

2 9.0180 20.0450 3.3639 

5 8.0748 17.5791 3.1646 

9.8281 21.9967 3.5091 

1 9.3763 20.7595 3.4303 
20 

2 8.9816 19.7049 3.3561 
5 8.0421 17.2877 3.1577 

9.7075 20.9726 3.4884 

1 9.2612 19.8083 3.4107 
10 

2 8.8713 18.8124 3.3375 
5 7.9434 16.5200 3.1415 



21 



4-3. Free vibration of plates 

Consider a plate of uniform thickness, h and with length and width as a 
and b, respectively. In this case, the displacement field is approximated 
by Lagrange elements (Q4 and Q8) and with NURBS basis functions. The 
computed frequencies for a square simply supported plate is given in Table |5j 
It can be seen that the numerical results from the present study are found to 
be in good agreement with the existing solutions. Table Ogives the results of 
non-dimensional frequency for first-order nonlocal plate theory for different 
values of nonlocal parameter, the plate aspect ratio and for various basis 
functions. Again, it can be seen that the nonlocal theory predicts smaller 
values of natural frequency than the local theory. 



Table 5: Comparison of fundamental frequency why/ p/G for simply supported plate. 



tt 



25j . mesh size 40 x40, *mesh size 8x8, 'Order of curve = 3, 5 x 5 control points. 



a/b a/h /i Ref. tf Method 



Q4 D Q8* NURBS* 



0.0930 0.0927 0.0926 0.0929 
10 1 0.0850 0.0847 0.0846 0.0849 

5 0.0660 0.0657 0.0657 0.0659 

1 

0.0239 0.0240 0.0238 0.0239 

20 1 0.0218 0.0219 0.0218 0.0219 

5 0.0169 0.0170 0.0169 0.0170 

0.0589 0.0588 0.0587 0.0590 

10 1 0.0556 0.0554 0.0554 0.0556 

5 0.0463 0.4620 0.0462 0.0464 

2 

0.0150 0.0150 0.0150 0.0151 

20 1 0.0141 0.0141 0.0141 0.0141 

5 0.0118 0.0118 0.0118 0.0118 

As another example, the nonlocal frequency -to- local frequency ratio (Qn/^l) 

22 



is computed for a simply supported isotropic square plate and the results are 
compared with those available in the literature 38|, |24( . Frequency ratio for 
different values of nonlocal parameter are presented in Table [6j It can be 
seen that the results from the present formulation are in good agreement 
with the existing results. 



Table 6: Comparison of frequency ratio for a simply supported square plate (a = 10, h 
0.34). tt 



25[, °[24|, D mesh size 40 x40, *mesh size 8x8, 'Order of curve = 3, 5 x 5 



control points. 



/i Ref.° Rcf.° 



Method 



Q4 



a 



Q8* 



NURBS* 



1.0000 1.0000 1.0000 1.0000 1.0000 

1 0.9139 0.9139 0.9099 0.9107 0.9107 

2 0.8467 0.8467 0.8393 0.8468 0.8468 

3 0.7925 0.7925 0.7857 0.7928 0.7928 

The influence of plate dimension on the frequency ratio for various nonlocal 
parameter is plotted in Figure OH]). The values of nonlocal parameter are as- 
sumed to vary between fi = (local elasticity) to fi = 3 nm 2 . As the length 
of the plate increases, the frequency ratio tend to increase and approach 
the local elasticity solution for considerably larger plate length, irrespective 
of the nonlocal parameter. The influence of nonlocal parameter is signifi- 
cant for smaller plate dimensions. Figure (Q shows the influence of nonlocal 
parameter on the frequency ratio for various plate dimensions. It can be 
seen that the nonlocal parameter has a stronger influence for smaller plate 
dimensions, as would be expected. To study the influence of nonlocal pa- 
rameter and the boundary condition on the natural frequency of an isotropic 
plate, a square plate of length 10 nm is considered with a/h = 100. The 
frequency ratio corresponding to two different boundary conditions, viz., all 
edges simply supported and all edges clamped boundary conditions are plot- 
ted in Figure (fTUj) . It can be seen that both the boundary condition and the 
nonlocal parameter has an influence on the frequency ratio. The higher the 
nonlocal parameter, the larger is the influence, irrespective of the boundary 
condition. 
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Figure 8: Effect of plate dimensions on the frequency ratio (VLnl/Ql) for a simply sup- 
ported square plate for various internal length with a/h = 100. 
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Figure 9: Effect of internal length on the frequency ratio (£Inl/Ql) for a simply supported 
square plate for various plate dimensions with a/h = 100. 
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Figure 10: Effect of internal length on the frequency ratio (Unl/^l) for all edges simply 
supported and all edges clamped boundary condition with a/h — 100 and a = 5. 



5. Conclusion 

In this paper, the constitutive model proposed by Eringen is used to model 
the Timoshenko beam and the plates based on FSDT. The natural frequencies 
of Timoshenko beam and first order shear deformable plate are studied using 
Lagrange polynomials and NURBS basis functions. Numerical experiments 
have been conducted to bring out the effect of boundary condition and the 
nonlocal parameter on the natural frequency of the beams and plates. The 
computed results are found to be in good agreement with analytical results. 
It can be seen that the NURBS basis functions requires fewer degrees of 
freedom to yield same order of accuracy as that of the RBF in case of the 
beams. It can be inferred that the effect of nonlocal parameter is to reduce 
the natural frequency of the plate irrespective of the boundary conditions. 
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